Biased bilayer graphene: semiconductor with a gap tunable by electric field effect 
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We demonstrate that the electronic gap of a graphene bilayer can be controlled externally by 
applying a gate bias. From the magneto-transport data (Shubnikov-de Haas measurements of the 
cyclotron mass), and using a tight binding model, we extract the value of the gap as a function of 
the electronic density. We show that the gap can be changed from zero to mid-infrared energies by 
using fields of ^1 V/nm, below the electric breakdown of Si02. The opening of a gap is clearly seen 
in the quantum Hall regime. 

PACS numbers: 81.05.Uw, 73.20.-r, 73.20.At, 73.21.Ac, 73.23.-b, 73.43.-f 



The electronic structure of materials is given by their 
chemical composition and specific arrangements of atoms 
in a crystal lattice and, accordingly, can be changed 
only slightly by external factors such as temperature 
or high pressure. In this Letter we show, both exper- 
imentally and theoretically, that the band structure of 
bilayer graphene can be controlled by an applied elec- 
tric field so that the electronic gap between the valence 
and conduction bands can be tuned between zero and 
mid-infrared energies. This makes bilayer graphene the 
only known semiconductor with a tunable energy gap 
and may open the way for developing photodetectors and 
lasers tunable by the electric field effect. The develop- 
ment of a graphene-based tunable semiconductor being 
reported here, as well as the discovery of anomalous inte- 
ger quantum Hall effects (QHE) in sing le layer P,[3l and 
unbiased bilayer Q graphene, which are associated with 
massless [3l and massive Q Dirac fermions, respectively, 
demonstrate the potential of these systems for carbon- 
based electronics [6|. Furthermore, the deep connection 
between the electronic properties of graphene and certain 
theories in particle physics makes graphene a test bed for 
many ideas in basic science. 

Below we report the experimental realization of a 
tunable-gap graphene bilayer and provide its theoretical 
description in terms of a tight-bindingjnodel corrected by 
charging effects (Hartree approach) [7[. Our main find- 
ings arc as follows: (i) in a magnetic field, a pronounced 
plateau at zero Hall conductivity axy = is found for 
the biased bilayer, which is absent in the unbiased case 
and can only be understood as due to the opening of a 
sizable gap, Ag, between the valence and conductance 
bands; [ii) the cyclotron mass, mc, in the bilayer biased 
by chemical doping is an asymmetric function of carrier 
density, n, which provides a clear signature of a gap and 
allows its estimate; {iii} by comparing the observed be- 



havior with our tight-binding results, we show that the 
gap can be tuned to values larger than 0.2 eV; {iv) we 
have crosschecked our theory against angle-resolved pho- 
toemission spectroscopy (ARPES) data [8| and found ex- 
cellent agreement. 

The devices used in our experiments were made from 
bilayer graphene prepared by micromechanical cleavage 
of graphite on top of an oxidized silicon wafer (300 nm 
of Si02) [§]. By using electron-beam lithography, the 
graphene samples were then processed into Hall bar de- 
vices similar to those reported in Refs. To induce 
charge carriers, we applied a gate voltage Vg between the 
sample and the Si wafer, which resulted in carrier con- 
centrations rii — aVg due to the electric field effect. The 
coefficient a = 7.2 x 10^°cm~^/V is determined by the 
geometry of the resulting capacitor and is in agreement 
with the values of ni found experimentally In or- 
der to control independently the gap value and the Fermi 
level Ep, the devices could also be doped chemically by 
exposing them to NH3 that adsorbed on graphene and 
effectively acted as a top gate providing a fixed elec- 
tron density no The total bilayer density n is then 
n = Til + no relatively to half-filling. The electrical mea- 
surements were carried out by the standard lock-in tech- 
nique in magnetic fields up to 12T and at temperatures 
between 4 and 300K. 

We start by showing experimental evidence for the gap 
opening in bilayer graphene. Figure [ija) shows the mea- 
sured Hall conductivity of bilayer graphene, which allows 
a comparison of the QHE behavior in the biased and un- 
biased cases. Here the curve labeled "pristine" shows 
the anomalous QHE that is characteristic of the unbi- 
ased bilayer 0]. In this case, the Hall conductivity ex- 
hibits a sequence of plateaus at Cxy = ANe^ /h where N 
is integer and the factor 4 takes into account graphene's 
quadruple degeneracy. The = plateau is strikingly 
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FIG. 1: (color online) (a) Measured Hall conductivity of 
pristine (undoped) and chemically doped bilayer graphene 
(no ~ 5.4 X 10^^ cm^^), showing a comparison of the QHE 
in both systems, (b) Cyclotron mass vs n, normalized to the 
free electron mass, mo. Experimental data are shown as O. 
The solid line is the result of the self-consistent procedure 
and the dashed line corresponds to the unscreened case. The 
inset shows an electron micrograph (in false color) of one of 
our Hall bar devices with a graphene ribbon width of 1 ^m. 



absent, so that a double step of 8e^/ft, in height occurs 
at n = 0, indicating a metallic state at the neutrality 
point . Note that the back-gate voltage induces asym- 
metry between the two layers but QHE measurements 
can only probe states close to Ep and are not sensitive 
to the presence (or absence) of a gap below the Fermi sea. 
To probe the gap that is expected to open at finite Vg, 
we first biased the bilayer devices chemically and then 
swept Vg through the neutrality point, in which case 
passes between the valence and conduction bands at high 
Vg. The energy gap is revealed by the appearance of the 

= plateau at axy = [see the curve labeled "doped" 
in Fig. [Ija)]. The emerged plateau was accompanied by 
a huge peak in longitudinal resistivity pxx , indicating an 
insulating state (in the biased device, pxx at n = ex- 
ceeded 150 kOhm at 4 K, as compared to w 6 kOhm for 
the unbiased case under the same conditions). The re- 
covered sequence of equidistant plateaus represents the 
"standard" integer QHE that would be expected for an 
ambipolar semiconductor with an energy gap exceeding 
the cyclotron energy. The latter is estimated to be > 40 
meV in the case of Fig. [Ija). 

To gain further information about the observed gap, 
we measured the cyclotron mass of charge carriers and its 
dependence on n. To this end, we followed the same time- 
consuming procedure as described in detail in Rcf. [ij 
for the case of single-layer graphene. In brief, for many 
different gate voltages, we measured the temperature 
(T) dependence of Shubnikov-de Haas oscillations and 
then fitted their amplitude by the standard expression 
T/ smh{2Tr^kBTmc/heB). To access electronic properties 
of both electrons and holes in the same chemically biased 
device, we chose to dope it to no ~ 1.8 x 10^^ cm~^, i.e. 
less than in the case of Fig. 1. The results arc shown in 



Fig.ffljb). The linear increase of rric with \n\ and the pro- 
nounced asymmetry between hole- and electron-doping of 
the biased bilayer are clearly seen here. 

To explain the observed Hall conductivity and cy- 
clotron mass data for bilayer graphene, in what follows, 
we shall use a tight-binding description of electrons in 
bilayer graphene. Its carbon atoms are arranged in two 
honeycomb lattices labeled 1 and 2 and stacked accord- 
ing to the Bernal order (^l-i32), where A and B re- 
fer to each sublattice within each honeycomb layer, as 
shown in Fig. [5] (a). The system is parameterized by 
a tight-binding model where TT-clcctrons are allowed to 
hop between nearest-neighbor sites, with in-plane hop- 
ping t and inter-plane hopping t±. Throughout the Let- 
ter wc use t = 3.1 eV and t± = 0.22 eV. The value of 
t is inferred from the Fermi-Dirac velocity in graphene, 
VP = {V2/i)at/h « 10^ms-\ where a « 2.46 A is the 
same-sublattice carbon-carbon distance, and t±_ is ex- 
tracted by fitting mc (see below). For the biased system 
the two layers gain different electrostatic potentials, and 
the corresponding energy difference is given by eV. The 
presence of a perpendicular magnetic field B = _B Cz is 
accounted for through the standard Pcicrls substitution, 
t ^ t cxp{ie J^^^ A • dr}, where e is the electron charge, 
5 the vector connecting nearest-neighbor sites, and A the 
vector potential (in units such that c = 1 = h). 

Fig. [2] (b) shows the electronic structure of the biased 
bilayer near the Dirac points {K or K'). In agreement 
with the Hall conductivity results in Fig. [TJa) , one can 
see that the unbiased gapless semiconductor (dashed line) 
becomes, with the application of an electrostatic poten- 
tial V , a small-gap semiconductor (solid line) whose gap 
is given by: Ag = [e^VHl/{tl + e^V^)]^^^. As V can 
be externally controlled, this model predicts that biased 
bilayer graphene should be a tunable-gap semiconduc- 
tor, in agreement with results obtained previously using 
a continuum model • Note that the gap does not reach 
a minimum at the K point due to the "mexican-hat" 
dispersion at low energies llj . 

The electric field induced between the two layers can 
be considered as a result of the effect of charged surfaces 
placed above and below bilayer graphene . Below is an 
accumulation or depletion layer inthc Si wafer, which has 
charge density ??ie. Dopants above the bilayer effectively 
provide the second charged surface with density noe. As- 
suming equal charge —ne/2 in layers 1 and 2 of the bilayer 
we find an unscreened potential difference given by. 



V = i2-n/no)noed/{2eo).. 



(1) 



where Eq is the permittivity of free space, and d ^ 
0.34 nm is the interlayer distance. A more realistic de- 
scription should account for the charge redistribution 
due to the presence of the external electric field. For 
given V and n, we can estimate the induced charge im- 
balance between layers An{n, V) through the weight of 
the wave functions in each layer (Hartrec approach; also. 
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FIG. 2: (color online) (a) Lattice structure of bilayer graphene 
and parameters of our model (see text), (b) Band structure 
of bilayer graphene near the Dirac points for eV — 150 meV 
(solid line) and V = (dashed), (c) eV as a function of 
n: solid and dotted lines are the result of the self-consistent 
procedure (see text) for t± — 0.2 eV and t± = 0.4 eV, respec- 
tively; dashed line is the unscreened result; circles represent 
eV vs. n measured by ARPES [1]. (d) Band gap as a func- 
tion of n (bottom axis) and Vg (top): solid and dashed lines 
are for the screened and unscreened cases, respectively. The 
thin dashed-dotted line is a linear fit to the screened result at 
small biases. 

see 0). This charge imbalance is responsible for an in- 
ternal electric field that screens the external one, and a 
self-consistent procedure to determine the screened elec- 
trostatic difference requires 

F = (2 - n/na + An{n, V)/no) noed/{2eo) . (2) 

Zero potential difference and zero gap are expected at 
n = 2no in both unscreened and screened cases, as seen 
from Eqs. (P) and ^ and the fact that An(n, 0) = 0. 

In Fig.[2](c) our calculations using Eqs. ([T]) and ([2]) are 
compared with ARPES measurements of the V depen- 
dence on n in bilayer graphene by Ohta et al. 0] . In their 
experiment, n-type doping with w lOx 10^^ cm~^ was 
due to the SiC substrate and therefore fixed. The elec- 
tronic density no induced by the deposition of K atoms 
onto the vacuum side was then used to vary the total den- 
sity. A zero gap was found around n « 23 x 10^^ cm^^ 
from which value we expect n}^ « 11 x 10^^ cm^^, in 
agreement with the experiment. In order to compare the 
behavior of V with varying n we replace no in Eqs. ([1]) 
and dll) with hq = n — n' . The result for the unscreened 
case [Eq. ([T|)] - shown in Fig. [5] (c) as a dashed line - 
cannot describe the experimental data. The solid and 
dotted lines are the screened results obtained with the 
self-consistent procedure [Eq. ^] for t± = 0.2 eV and 



t± = 0.4 eV, respectively; both are in good agreement 
with the experiment, except in the gap saturation regime 
at n > 50 X 10i2cni-2. 

For the experiments described in the present work, the 
expected behavior of the gap with varying n or, equiv- 
alently, Vg is shown in Fig. [2] (d). The dashed line is 
the unscreened result [V given by Eq. ([J)] and the solid 
line is the screened one [Eq. ([2])]. In both cases, the 
chemical doping was set to no = 1.8 x 10^^ cm^^ at 
which rUc was measured in our experiment (equivalent of 
Vg « 25 V). The dashed-dotted (blue) hne is a linear fit 
to the screened result for small gap, yielding Ag(meV) ~ 
(3\Vg{Y) - 25| with a coefficient /3 w 1.2 meV/V. The lin- 
ear fit is valid in the small-gap regime (Ag ^ t±) only, 
and the theory predicts a gap saturation to Ag ^ t± 
at large biases. Note that the breakdown field for Si02 
is 1 V/nm (i.e. 300 V for the used oxide thickness) and, 
therefore, practically the whole range of allowed gaps (up 
to t±) should be achievable for the demonstrated devices. 

To explain the observed behavior of the cyclotron 
mass, rrif., shown in Fig.[T]^b), we used the semi-classical 
expression mc{n) = {fi^ /2Tr)dA{E)/dE\E=EF{n)^ where 
A{E) is the fc-space area enclosed by the orbit of energy 
E and n the carrier density at Ep. In Fig.[ljb) our the- 
ory results are shown as dashed and solid lines for the 
unscreened and screened description of the gap, respec- 
tively (analytical expressions for rUc in the biased bilayer 
will be given elsewhere [l^ ) . The inter-layer coupling t± 
is the only adjustable parameter, as t is fixed and V is 
given by Eq. ([T]) or Eq. ([2]) . The value of t± could then 
be chosen so that theory and experiment gave the same 
TOc for n « 3.6 x lO^^cm"^. At this particular density the 
gap closes and the theoretical value becomes independent 
of the screening assumptions. We found t± ~ 0.22 cV, in 
good agreement with values found in the literature. The 
theoretical dependence mc(n) agrees well with the exper- 
imental data for the case of electron doping. Also, as seen 
in Fig.[Ub), the screened result provides a somewhat bet- 
ter fit than the unscreened model, especially at low elec- 
tron densities. This fact, along with the good agreement 
found for the potential difference data of Ref. Q [see 
Fig. [2] (c)] , allows us to conclude that for doping of the 
same sign from both sides of bilayer graphene, the gap 
is well described by the screened approach. In the hole 
doping region in Fig. [IJb), the Hartree approach under- 
estimates the value of TOc whereas the simple unscreened 
result overestimates it. This can be attributed to the fact 
that the Hartree theory used here is reliable only if the 
gap is small compared to t±. In our experimental case, 
no > and, therefore, the theory works well for a wide 
range of electron doping 7i > 0, whereas even a modest 
overall hole doping n < corresponds to a significant 
electrostatic difference between the two graphene layers. 
In this case, the unscreened theory overestimates the gap 
whereas the Hartree calculation underestimates it. How- 
ever, it is clear that the experimental data in Fig. [Ub) 
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FIG. 3: (color online) Energy spectrum for a ribbon of bilayer 
graphene with zigzag edges, t^/t — 0.2, B = SOT, and width 
TV = 400 unit cells: (a) eV = 0; (b) eV = tx/10. (c) Sketch 
of the bands close to zero energy (for the biased bilayer) with 
indication of bulk (solid lines) or edge (dotted lines) states 
and their left (L) or right (R) positions along the ribbon. 
Quasidegeneracies have been removed for clarity. 



interpolate between the screened result at low hole dop- 
ing and the unscreened one for high hole densities. This 
indicates that the true gap actually lies between the un- 
screened and screened limits [see Fig. [2](d)], and that a 
more accurate treatment of screening is needed when eV 
becomes of the order of t± . 

In what follows, we model and discuss the QHE data 
presented in Fig. [TJa). We consider a ribbon of bilayer 
graphene [ist with zigzag edges (armchair edges give sim- 
ilar results). Fig. [3] shows the energy spectrum in the 
presence of a strong magnetic field. Panel (a) corresponds 
to the unbiased case [see the curve labeled "pristine", 
Fig.IUa)], where the four degenerate bands at zero energy 
contain four degenerate bulk Landau levels Q and four 
surface states characteristic of the bilayer with zigzag 
edges [lH • The spectrum for a biased device is shown in 
panel (b). In this case two flat bands with energies —eV/2 
and eV/2 appear, similar to the case of zero magnetic 
field. The other two zero energy bands become dispersive 
inside the gap, showing the band-crossing phenomenon. 
The Landau level spacing is set by 7 = ^/ijltajl b (^b is 
the magnetic length), and as long as eV <C t±, the bias is 
much smaller than the Landau level spacing at low fields. 
Then non-zero Landau levels in the bulk are almost in- 
sensitive to V, as seen in Fig. [3] (b), except for a small 
asymmetry between Dirac points. A close inspection of 
Fig. [3] shows that the valley degeneracy is lifted due to 
the different nature of the Landau states at K and K' 
valleys with respect to their projection in each layer. The 
valley asymmetry has a stronger effect in the zero energy 
Landau levels, where the charge imbalance is saturated. 
This opens a gap of eV in size. Also, there is an intra- 



valley degeneracy lifting [Fig.[3](c)], because only one of 
the two Landau states of the unbiased system remains 
an eigenstate when a bias is applied. For eV > t±_ (not 
shown in Fig. [3]) the dispersive modes start crossing with 
non-zero bulk Landau levels. 

Let us now model the measured Hall conductivity for 
the biased bilayer graphene, which is shown in the (red) 
curve labeled "doped" in Fig. [2a) . We consider the case 
of the chemical potential lying inside the gap, between 
the last hole- and the first electron- like bulk Landau 
levels, and crossing the dispersive bands as shown in 
Fig. [3](c). As pointed out by Laughlin changing the 
magnetic fiux through the ribbon loop by a flux quan- 
tum makes the states to move rigidly towards one edge. 
In the usual integer QHE, the energy increase due to 
this adiabatic flux variation results in the net transfer 
of n X g electrons (spin and valley degeneracy g) from 
one edge to the other, and the quantization of the Hall 
conductivity follows the expression 1/V = gne^ /h, 
where X is the current carried around the loop and V 
the potential drop between the two edges. However, in 
the present case there is no net charge transfer across 
the ribbon. As seen in Fig. [3] (c), the band states at 
the chemical potential belonging to the same band are 
surface states localized at the same edge (see the figure 
caption for details). The rigid movement of the states to- 
wards one edge makes an electron-hole pair to appear at 
both edges, resulting in zero net charge transfer. There- 
fore, we expect a Hall plateau Uxy = showing up when 
the carriers change sign, i.e. at the neutrality point. Ac- 
cordingly, the Hall conductivity of the biased bilayer is 
given by axy = 4:Ne'^/h for all integer N, including zero. 
Note that at the zero Hall plateau the current carried 
around the ribbon loop is zero, X = 0, which implies, 
from the theory view point, a diverging longitudinal re- 
sistivity at low T, in stark contrast to all the other Hall 
plateaus that exhibit zero pxx-, as in the standard QHE. 
This behavior has been observed experimentally, as dis- 
cussed above with reference to Fig. [Ija). This concludes 
our interpretation of the experimental data. 
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